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I. Introduction 

The Lee-Friedrichs model 1 ' 2 has been very useful for the study of the properties of 
unstable systems. It is completely soluble, and provides a framework for the analytic study 
of decay and scattering systems 3 . The construction of generalized states with exact expo- 
nential decay 4 in the rigged Hilbert space 5 (Gel'fand triple) was originally motivated by 
the properties of this model; such states have found application in the theory of iterated 
maps 6 and play an important role in the study of irreversible processes 7 . For both theoret- 
ical development and physical applications, it is important to have a relativistic model for 
unstable systems. An unstable particle generally decays into a final state of two or more 
particles in a process for which the total energy is conserved, but the total mass is not. 
The equivalence between mass and energy, which enters quantitatively in the kinematical 
description of such a process, is a fundamentally relativistic relation. We shall use a for- 
mulation of relativistic quantum theory (and quantum field theory, as will be discussed in 
Section 2) which is well-suited for this problem in that it treats the energy and momentum 
of a particle as independent observables; hence, the mass becomes a dynamical variable. 
The transition from the undecayed initial state to the decayed system in the final state 
can then be studied as a continuous quantum mechanical evolution. 

This formulation 8 , a generalization to the N— body system of the relativistic quantum 
theory of Stueckelberg 9 developed for a single particle, describes evolution with a Poincare 
invariant parameter r, sometimes called the "historical time." This invariant parameter 
corresponds to the universal time of (freely falling) ideal clocks 10 , and may be identified 
with the universal time originally postulated by Newton. In this quantum theory, the 
four momenta p^ are conjugate to the four coordinates of space and time x^ = (t,x). 
Classically, t, the time at which an event occurs in the laboratory (measured on a standard 
ideal clock, i.e., running with r), is a dynamical variable on the same footing as x (the place 
at which the event occurs in the laboratory). The wave function, for example, for a one- 
particle system, ip T (x) (we shall use x to represent x^), is the amplitude for the probability 
density |^ T (x)| 2 to find an event at (t, x) at the historical time r. Since this density is 
normalized with respect to integration over t, x (i.e., it belongs to the Hilbert space L 2 (R 4 )) 
measurable unbounded excursions in t, as for x in the non-relativistic quantum theory, are 
not admissible 11 . 

Classically, the variables t(r), x(r), correspond to the value of t, x of a signal detected 
in the laboratory (in some local inertial Lorentz frame) for an event generated at time r 
(in some, perhaps other, frame). The result of the detection experiment is influenced by 
forces through the Hamilton equations (i = 1, 2, . . . , N for an N— body system) 



where K is a function on the 8iV dimensional phase space, as well as the state of motion 
of the frame by means of the Lorentz transformation. 

In the quantum theory, the Stueckelberg-Schrodinger equation for an N— body system 
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dp^ 
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(1.1) 
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(X 1 . ..X N ) = K(x 1 . ..X N ,p 1 . ..p N )^ T (x 1 . ..X N ) 



(1.2) 
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describes the evolution of the wave function, and the dynamical variables satisfy the Heisen- 
berg equations 

df=^ = i[K,x'>] p» = ^ r = i[K,p% (1.3) 
with [g^ v = (— , +, +, +); we shall use units for which h = c = 1] 

[ X »,p»]=ig^. (1.4) 

For a single free particle, the assumption 

K = (1.5) 

where M is an intrinsic parameter belonging to the particle which sets the scale of energy- 
momentum relative to its space-time evolution in r, results in (both classically and quan- 
tum mechanically) 

~d~r = ~M' (L6) 
For the classical case, one may eliminate dr to obtain 

— = - (17) 
dt E' [ } 

the Einstein relation for velocity. A localized ( in x, t) free wave packet moves with r, 
according to the Ehrenfest theorem, approximately along the classical world line. The 
generic form (one may also include electromagnetic interactions in gauge invariant form) 

N p 

^ = E^£r+ v( Xl ...z N ), (i.8) 

where V is a Poincare invariant function of the {xf }, provides a family of consistent di- 
rect action potential models to describe the evolution of a relativistic iV-body system. 
The Gibbs ensembles for the equilibrium relativistic statistical mechanics of such systems, 
both for distinguishable and indistinguishable particles has been worked out 12 , as well 
as the BBGKY and Boltzmann equations 13 for the non-equilibrium evolution of the iV- 
body system. The equilibrium limit of the Boltzmann distribution provides a relativistic 
equilibrium mass distribution, for which the low energy and non-relativistic limits have 
been investigated (with results, in this limit, of course, in agreement with non-relativistic 
statistical mechanics), and the relativistic Fermi and Bose equilibrium distributions are 
currently being studied, for example, for their consequences in astrophysics and condensed 
matter theory 14 ' 15 . Van Hove 16 has used the non-relativistic Lee-Friedrichs model to dis- 
cuss some of the fundamental properties of non-equilibrium statistical mechanics, and we 
shall, in a future publication, apply the model that we have developed here to analogous 
considerations. 
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As an illustration of the applicability of the model (1.8), the two-body problem with 
"central" potential has been studied in detail 17 . Since we shall make use of two body 
kinematics in the study of the decay problem, we shall review here some of the details of 
the treatment of this problem. For the two-body case, we take 

where p = ^x^x^, and the relative space-time coordinates are given obey 

x M = a ,M_ a ,M ( 110 ) 

In the non-relativistic limit, t\ — > t^-, so that p — > r = |xi — X2I. The total momentum of 
the system is 

P»=p»+p%, (l.H) 

which commutes with the relative coordinate x^. The operator K, as for the non- 
relativistic analog, can therefore be written as 

pfi p 

K = + Krel, (1.12) 



where 

K r 

„ _ Mtf - M lP % (L13) 
P M 1 +M 2 ' 

corresponds to the evolution operator for the relative motion. Clearly, K has only contin- 
uous spectrum, since the center of momentum motion is that of a free particle. Defining 

x* = M^1±M^1, (1 . 14) 

where 

M = Mi+M 2 , (1.15) 

we see that by (1.1), 

— = —. (1.16) 

However, K re i, of (1.13), can have discrete or continuous spectrum in spite of the fact 
that p^p^ is a hyperbolic operator. One studies this problem by separating the center of 
momentum variable from K (i.e., taking a direct integral representation of the solution 
over P M , the absolutely conserved total energy-momentum). It is a remarkable fact that 
if we take the support of the eigenfunction V , fc( a; )? °f the equation 

K re iip k (x) = kil) k {x) (1.17) 
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to lie in the sub-measure space defined by 

x\ + x\ - t 2 > 0, (1.18) 

called the RMS (restricted Minkowski space), the spectrum of K re i then coincides with 
the spectrum of the non-relativistic Schrodinger equation for all central potentials V(r) = 
V(p)\ p=r . The corresponding energies are then determined by (1.12), i.e., in the center of 
momentum frame, for which P = 0, 

E = y/2M(k-K). (1.19) 



With the assumption that the particles have masses close to Mi, when they are asymp- 
totically spacelike separated (for which V ~ 0), and that the corresponding evaluation of 
K is constant over the entire spectrum of K re i [we argue that t-dependent perturbations 
of V at long intervals can shift the spectrum of K re i, but leave the absolutely conserved 
K constant], (1.19) becomes 

E = y / 2Mk + M 2 . (1.20) 

For the case of discrete spectrum, where k = k n , n integer, when the excitations are small 
in magnitude compared to M, 

En ^M + k n + MO((^) 2 ). (1.21) 

The same decompostion can be carried out for continuum values of k (scattering states) 
small compared to M. 

The RMS has a natural coordinatization in terms of the compact parameters 9, cf> and 
the non-compact parameter (3 (along with p) as 

x\ = p cos 4> sin 9 cosh Q 

xo = p sin d> sin 9 cosh 3 

n (1-22) 
x = t = p sin 9 sinh 8 

X3 = pcos6>; 



clearly, x\ + x\ - t 2 = p 2 sin 2 9 > 0, and x\ + x\ + x\ - t 2 = p 2 > as well. This 
submeasure space is 0(2, 1) invariant, but not 0(3, 1) invariant. The representations of 
0(3, 1) are induced from the irreducible representations of 0(2, 1) which are provided by 
the eigenfunctions of (1.17) 17 . 

A model for relativistic hydrogen is obtained by taking V = —e 2 /p. The exact eigen- 
functions are therefore known for this problem; it has recently been shown that the selec- 
tion rules for radiation coincide with those predicted by the non-relativistic theory 18 , and 
a relativistic form of the Zeeman effect has been worked out 19 . Many interesting group 
theoretical problems have been posed by this structure. The dynamical group for the rela- 
tivistic Kepler problem has been found 20 , and a general investigation of such non-compact 
dynamical groups is being carried out 21 . 
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The theory of unstable systems that forms the basis of our discussion in this work was 
first given in a systematic way by Wigner and Weisskopf 22 , in their work on the spectral 
line width for atomic radiation. In this (non-relativistic) theory, the unstable system is 
represented by a wave function ip which is usually taken to be an eigenfunction of an 
unperturbed Hamiltonian Hq; the action of the full perturbed Hamiltonian H then causes 
a non-trivial evolution of the wave function ip t = e~ lHtr ip, and the amplitude 

A(t) = (i/>,e" iH V) (1-23) 

is understood to be the amplitude for survival of the initial state, i.e., 

p{t) = \A{t)\ 2 (1.24) 

is the probability to find the initial state "undecayed" after a time t. After a time t, not too 
short and not too long, the evolution (1.23) is, for a wide class of models, well-approximated 
by the exponential of a complex number, i.e., of the form exp{— i(E — i^)t}, where 7 is 
the exponential decay rate for p(t). If H is defined on ip, it is well known that the survival 
probability p(t) has zero derivative at t = 0, and, moreover, it goes asymptotically, for 
large t, as an inverse power 23 of t. 

It should be remarked in this connection that a fundamental assumption underlying 
the Wigner- Weisskopf theory is that, under continuous evolution, 

e~ im ^ = A(t)1> + , 

where ip 1 - is the part of e~ tHt/ ijj orthogonal to ip. The linear superposition A^ip + ip 1 - is, in 
fact, not seen in the laboratory. One generally sees only the state ip or the decay products. 
It appears that the time of decay is a superselection rule, and that this linear superposition 
does not correspond to the physical situation. This problem has been discussed by Horwitz 
and Piron 24 , where an alternative formulation of the problem for the non-relativistic case 
has been given. We shall not discuss the relativistic generalization of that treatment here, 
but confine ourselves to the usual Wigner- Weisskopf form which, although perhaps not the 
most general theoretical framework, has led to a good quantitative approximation for the 
non-relativistic theory. 

One may understand the time evolution predicted by the Wigner- Weisskopf theory 
most easily by examining the inverse Laplace transform 25 (one takes A(t) = for t < 0) 

A(t) = ±-[ dz(^, - J_V0e-*** (1-25) 

Z7TZ J (j Z — tl 

where C is a contour that depends on the nature of the spectrum of H. 

If H has continuous spectrum, for example, on (0, 00), then the function 

R(z) = (V, -^-=V0 (1-26) 
z — tl 

is analytic in the upper half plane along with its extension in the first Riemann sheet to 
the lower half plane, excluding the cut along the positive real axis. The contour C then 
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corresponds to a line running above the positive real axis from oo to 0, then going around 
the branch point and continuing from to oo below the real axis. The lower part of this 
contour can be deformed to a vertical line from to —zoo for t > 0, providing only a small 
contribution for large t. The part of the contour lying above the real axis may then be 
deformed through the cut to the second Riemann sheet, lowering it eventually to a line in 
the second sheet from —zoo to 0, going around the branch point to connect with the part 
of contour that remains in the first sheet. In the course of this deformation, the contour 
passes the pole which has been shifted (for small coupling) from the discrete (bound) initial 
state embedded in the continuum on the real axis into the second sheet in the lower half 
plane, and this pole and its residue, proportional to e~ l< * t , where ( is the position of the 
pole, corresponds to the exponential decay law. For large t, this contribution becomes 
small compared to that provided by the neighborhood of the branch point, which yields an 
inverse power of t. For very short times, the contribution of the deformed integral cannot 
be ignored, and, in fact, the integral is best approximated along the real axis 26 . Since H 
is self-adjoint, clearly dp(t)/dt vanishes at t = 0, and it is easily seen that the power series 
in t starts as 1 — 0(t 2 ). 

Friedrichs 2 proposed the following soluble model, which provides a closed analytic 
form for what is called the reduced resolvent R(z), for which the analysis described above 
may be carried out in a very simply. Let 

H = H + V, (1.27) 

where Hq has an absolutely continuous spectrum in (0, oo) and a single bound (discrete) 
state ip for which 

H ij> = E^. (1.28) 
Furthermore, suppose that on the continuous spectrum of Hq characterized by 

(E\H \f) = E{E\f), (1.29) 

for any / in the domain of Hq, and that 

(E\V\E') = (1.30) 

for all E,E'. The diagonal part (ift\V\ip) of V, leading to a simple shift in Eq, is usually 
taken to be zero. Then, it follows from the identity (second resolvent equation) 

1 1 (I*) 



z — H z — Hq z — Hq z — H 
that 

R(z) = -J— + -J— I dE ^\V\E)(E\^—\^) (1.32) 

Z — Eq Z — Eq J Z — ti 

and 

(SI-J—IVO = -^=(E\V1,)R{z). (1.33) 
z — ti z — E 
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The simplifying truncation in (1.33) follows from the essential feature (1.30) of the 
model, implying physically that there are no "final state" interactions. This is often, in 
fact, a good physical approximation in decaying systems. Substituting (1.33) into (1.32), 
one finds 

h(z)R(z) = 1, (1.34) 



where 



dE 



\{E\V\1>) 



E 



h(z) = z — Eq - 

The singularities of R(z) arise from the vanishing of the function h(z). Since 



(1.35) 



ImH(z) = Imz 



1 + dE 



\z-E\ 



(1.36) 



there can be no zero in the first sheet of the function for Imz ^ 0. For sufficiently small 
V, it is easy to show that there is no zero on the negative real axis. Passing to the second 
sheet continuously through the cut, one finds 25 



h n (z) = h(z) + 2niW(z), 
where W(z) is the analytic extension of 

W(E) = \(E\V\^ 12 



(1.37) 



to a region of the lower half plane sufficiently large [that is, we assume an analytic function 
in this domain whose boundary is W(-E)] to include the possible zero of h n (z). If this 
zero occurs for z = C with small imaginary part, then 

= -27r|(£ |W)| 2 , (1.38) 

so that, recalling the representation 

lim- e = S(x), 

one sees that the pole occurs at 

ImC = -7r|(£ |W)| 2 . (1.39) 



1+ dE 



\{E\vm 

\z-E\ 2 



It is clear that for the description of, for example, particle decay, involving a real 
change in the the total particle mass of the system, a relativistic theory is required. Since, 
as pointed out by Henley and Thirring 27 , in the standard on-shell relativistic quantum 
field theory, both the annihilation operator for a particle and the creation operator for an 
antiparticle occur linearly combined in the local field operator, a relativistic analog of the 
model described above (in the framework of quantum field theory, in the form given by 
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Lee 1 ) would contain, along with the process n <-> p+ir~ , automatically and simultaneously, 
the process p <-> n + 7r+. Zumino 28 has remarked that the expression E = a/p 2 + to 2 is not 
analytic, and creates singularities that make the analytic continuations we have discussed 
above difficult; he suggested, toward the end of his article, the use of the quadratic form 
p^Pft instead, a form which we will, in fact, utilize in the formulation we give in the next 
sections. In the field theoretic form, the "coupling" which plays the role of (E\V\ip), 
is a function of the momentum of one of the particles emitted in the decay; its Fourier 
transform to configuration space results in a non-local function in space which would make 
Lorentz invariance difficult to achieve. The original form of the model proposed by Lee 1 , 
which was constructed to achieve the simplest form, did not allow for recoil, and hence its 
relativistic generalization would destroy Poincare invariance; however, in the form given 
by Van Hove 29 , recoil is taken into account. It is for this form of the field theory which we 
shall provide a relativistic generalization. 

It is the purpose of this paper to show how such a relativistically covariant version 
of this model can be constructed. In Section 2, we shall formulate the decay problem 
in terms of the relativistic quantum field theory resulting from the second quantization 
of the covariant quantum theory we have discussed above. It has recently been shown 
that such a field theory results in a consistent quantum electrodynamics 30 involving a five 
dimensional generalization of the Maxwell field. This field (including a Lorentz scalar 
field as the fifth component) reduces to the usual Maxwell field when a characteristic 
correlation length reaches a critical value. The theory has been shown 31 to provide insight 
into regularization (such as that of Pauli-Villars 32 ' 33 ) and renormalization methods used in 
standard quantum electrodynamics (e.g., refs. 34). We shall not describe here this general 
theory, but use the relativistic matter fields of this theory in much the same way that Lee 1 
and, in particular, Van Hove 29 used the non-relativistic quantized fields to describe the 
soluble Lee model. The information we obtain from this field theory, in a model which 
has two conserved quantum numbers defining particle number sectors, is then used, in 
Section 4, to construct a corresponding model in a single Hilbert space (which is not a 
tensor product Fock space representing a quantum field theory) according to the method 
of Friedrichs 2 . We thus achieve a relativistically covariant Lee-Friedrichs soluble model for 
an unstable system. 

It is a pleasure to dedicate this paper to Professor Jean-Pierre Vigier on the 45 th 
anniversary of his joining Louis de Broglie at the Institut Henri Poincare. The theory 
of E.C.G. Stueckelberg provides a framework in which the notion of de Broglie 35 , that 
the mass of a particle is associated with a quantum frequency, is realized through the 
identification of the energy in a relativistic theory as an observable kinematically indepen- 
dent of the momentum. The work of Vigier and collaborators (ref. 8), contributing to 
the development of this theory, we feel, makes the present work appropriate for such an 
occasion. 

2. Field Theory Model 

In order to determine the structure and physically motivated basic assumptions to be 
imposed on a Friedrichs type model for the relativistic description of an unstable system in 
a "one- particle" Hilbert space, we first study in some detail a model in relativistic quantum 
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field theory. We shall assume boson matter fields in what follows. The properties of the 
quantum fields may be obtained by a variety of methods; one of these is the construction of 
a Lagrangian density for which the variation of the associated action yields a field equation 
which coincides with the Stueckelberg-Schrodinger equation. A Lagrangian density which 
yields this result is, for example 30 , 



C = \{i^ T d T ^ T - d^ T d^ T + C.C.} 

— / dx dx' i(jI(x)i(jI(x')V(x — x')ilj T (x')ifj T (x)}, 



and the action is defined as 



= J dr J d 4 x C. (2.2) 



It then follows that the field conjugate to ip is IT = iip^ , so that the equal-r commu- 
tation relations are 

[^ T {x)^ T {x')] = 5\x-x'). (2.3) 

In momentum space, for which 

MP) = j^y J d 4 x e-^ x »M*)> ( 2 - 4 ) 

this relation becomes 

[Mp)^r(p')] = t5 4 (p-p'). (2.5) 

To see how this relation goes over to the on-mass-shell field theory, we follow a non- 
rigorous but instructive procedure. Multiply both sides of (2.5) by AE, and notice that 
for E = a/p 2 + m 2 , where m 2 is a dynamical variable, AE = Am 2 /2E, we obtain 

[Vv(p), flip')] = 2iE5 3 (p-p'), (2.6) 



where Vv(p) = y Am 2 i/) T (p), and n r (p) = y Am 2 U T (p). We may then take the limit 
Am 2 — > 0, m 2 — > mo, a fi xe d constant (in this limit, t — > r as well 8 ). The operator-valued 
distribution ip T (p) therefore peaks in support (in all matrix elements) sharply at m 2 ~ t/Iq 
in this limit. Eq. (2.6) is the conventional equal time commutation relation. 

The example we have studied above assumes positive energy. Both fields ^II may 
have negative energy as well, but the right hand side vanishes if one has positive, and 
the other negative (AE5(E — E') would be zero in this case). We remark, furthermore, 
that according to (2.5), the field ip(x) contains only the annihilation operator ip T (p), for 
both positive and negative energy, but not a creation operator for the antiparticle (as 
for the second quantization of the Schrodinger wave function). The CPT conjugate of a 
negative energy function corresponds to the representation of the antiparticle observed in 
the laboratory. This is consistent with Stueckelberg's (and Feynman's) 9 original interpre- 
tation, associating the antiparticle with the particle going "backwards" (with r) in time. 
The quantum electrodynamics (a U(l) gauge field in interaction with the matter field) 
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associated with this theory is obtained by requiring gauge invariance of the Stueckelberg 
equation 30 . In addition to the derivative —id 11 , which becomes the gauge covariant deriva- 
tive —id^—a^(x, r), the derivative id T also requires compensation, i.e., id T — > id T +as(x, r), 
introducing a fifth, Lorentz scalar, field. 

It has been shown that the Maxwell potential (satisfying the usual Maxwell equations 
with a four dimensional conserved current as its source) is obtained from the field a M (x, r) 
by integration over r 30 , i.e., the Maxwell field is the zero mode of the "pre-Maxwell" field 
(the a 5 field decouples in the U(l) theory). In the limit in which the Fourier transform 
a^{x, s) on r of the pre-Maxwell fields have support only in a small interval As near zero (a 
long correlation length limit) one can show that the theory goes over to the usual form of 
quantum electrodynamics 30 . The theory has been consistently quantized, using canonical 
methods 36 , in full generality. In the quantized form, the theory goes over to the standard 
quantum electrodynamics, carrying Pauli-Villars regularization 33 ' 34 , in the long correlation 
length limit, as we have remarked above. We shall not discuss this structure further here, 
but concentrate on a model field theory with a generalized invariant Hamiltonian of the 
form 

K = K + Ki (2.7) 
where (we write p 2 = p^p^ , k 2 = k^k^ henceforth) 

Ko = J d 4 p^-tf(p)b(p) 

+ J d 4 p^aUp)a N ( P ) (2-8) 

+ / d 4 k^-a^{k) e a e {k\ 

and 

Kj = ! d 4 k d 4 p {f{k)b\p)a N {p - k)a e (k) 

J (2.9) 

+ f(kr^ N (p-k)al(k)b(p)}, 

describing the process V <-> N+9. Here, b(p) is the annihilation operator for the V particle, 
and cin(p), a$(k), for the N and 9 respectively. This model is a relativistic generalization 
of Van Hove's form 29 of the Lee model for the non-relativistic theory (note that expressing 
the "kinetic energy" in the on-mass-shell form a/p 2 + rriQ would make it impossible to 
separate the center of momentum motion in the two-body channel). The operators 

Qi= f d 4 p [6t(p)6(p) + aUp)a N (p)} 

J (2.10) 
Q 2 = d 4 p [a j N (p)a N (p) - al(p)a e (p)] 



are conserved (as one easily sees using the commutation relations (2.5)), and enable us to 
decompose the Fock space to sectors. We shall study this problem in the lowest sector 
Qi = 1, Q 2 = for which there is just one V or one iV and one 9. 
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We now wish to calculate 

tf T = e~ iKT ^ , (2.11) 
where we take the normalized wave packet 

*o = J g(p)b\p)d*p\0) (2.12) 

as the one-particle initial state (the ^-particle). We can now compute the relativistic 
analog of the survival amplitude 

A a (r) = (* , *r) = (*o, e- iKT V ) (2.13) 

as a function of the invariant time r. The initial state may be chosen so that r is the mean 
proper time for the undecayed system (i.e., for which ((dx^ / dr)) {(dx M / dr)) = —1). 

Since A(r) is Lorentz invariant, we may transform from an arbitrary frame to the 
(average) rest frame of the initial state under the evolution induced by K. Its r-dependence 
therefore clearly can arise only from the evolution of to a different channel (assumed 
here to be more rapid than the significant spreading of the wave packet). 

In the sector Q\ = 1 and Q 2 = 0, conserved by the evolution, a general representation 
for the Fock space state \I/ r is 

¥ T = J d 4 pA(p,rp(p)\0} + J d 4 pd 4 k B(p,k,r)aUp)al(k)\0), (2.14) 

where we choose for the initial condition 

A( P ,0) = g( P ) 
B(p,fc,0) = 0, 

so that the initial state is that given by (2.12). The evolution of the system is given by 

i ^=KV T , (2.16) 

the quantum field theory form of Eq. (1.2). Applying the definitions (2.8), (2.9) for K 
and the commutation relations for the annihilation-creation operators, one obtains the 
differential equations 

ifrip, r) = l^MP, r)+J d 4 kf(k)B(p - k, k, t) (2.17) 



and 



3B f n 2 k 2 \ 

1 a7 fe *• T) = \ir N + 2M e ) Bip - k ' T) + r{k)A{ ? + k ' T) - (2 ' 18) 
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Defining the Laplace transforms 

e lZT A(p, r)dr 

r oc (2-19) 

B(p,k,z)= / e izr B( Pl k,r), 
Jo 

we obtain from (2.17), (2.18) the relations 

[z - A(p, z) = ig(p) + J d*kf(k)B(p - fc, fc, z), (2.20) 

where the term igf(p) arises from the boundary point at r = of the integration by parts 
(the contribution at r — > oo vanishes for 2; in the upper half plane), and 

*5(P, k, z) = + ^)B(p, k, z) + r(k)A(p + k,z). (2.21) 

Substituting (2.21) into (2.20), we obtain the covariant analog of Eq. (1.34), 

h(p,z)A(p,z)=ig(p), (2.22) 

where 

•to* J*" ; (2 ' 23) 

^ 2M N 2M g 

In this form, one recognizes the essential content of the relativistic Lee model. The 
term p 2 /2My corresponds to the mass (squared) of the unstable system, i.e., the discrete 
(for fixed p M ) eigenvalue of the unperturbed initial state, and the cut, the singularity 
corresponding to the continuous spectrum of the decay channel. 

The Laplace transform of the survival amplitude (2.13), 

A s (z) = J™ e lZT A s (r)dr = i(* , j^*o), (2.24) 
analytic for z in the upper half-plane, is then given, according to (2.12) and (2.14), by 

A s (z) = J d 4 pg(pyA(p,z), (2.25) 



i.e.(with (2.22)), 



The inverse transform is 



(,) = */ '^g. (2.26) 



A s (r) = ^-J G e-^A s (z)dz, (2.27) 
13 



where C is a contour taken from +00 + it to —00 + ie, i.e., just above the real axis. In the 
case of semi-bounded spectrum for the continuum channel, for r > 0, the left hand part of 
the contour can be deformed to the lower half-plane, where this part of the contour (except 
for the neighborhood of the branch point) provides a small background contribution for 
large r. In the same way as discussed in Section 1, the right hand half line can be rotated 
into the second sheet, and brought down to the line (0, 00), picking up the contribution of 
the pole of h(p, z) on the way. If the spectrum is unbounded, analytic continuation can be 
made by bringing the whole line of integration down into the lower half plane. 

Note that the contribution of the pole corresponding to the zero of the analytic con- 
tinuation of h(p, z) is spread out by the integration (2.26), but this spreading can be very 
small (it occurs primarily in the real part of the pole, since the imaginary part, for weak 
coupling, is almost independent of p) . In the next section, we discuss the structure of the 
spectrum of the continuous channel. 



3. The Spectrum 

Clearly, the unperturbed spectrum (for f(k) = 0) for the relativistic Lee model, 
according to (2.8), is purely continuous. We see, however, from the structure of (2.23) 
that for each value of p 2 , there is an effective discrete eigenstate with eigenvalue p 2 /2M V: 
considering the Hilbert space as a direct sum space over the absolutely conserved total 
energy-momentum, which coincides with of the system. The remaining terms in (2.8) 
can be written as 

/ ^^{^^ + ^} a lv(P- fc )4(^(^(P-^) 5 (3-1) 

since we can replace p — k by p; the integration of the first term over k then provides a 
factor which is the number operator of the 6 field, and the second term has a factor, after 
integration over p, which is the number operator for the N field. In the sector Qi = 1 
and Q2 = 0, these factors are unity in the continuum channel and zero in the initial state 
channel, and hence do not change the valuation of this contribution to the unperturbed 
spectrum. Clearly, (3.1) corresponds to the continuous spectrum associated with the cut 
in (2.23), for each fixed p M . We therefore study the range of the function 

(p-k) 2 k 2 , N 

The range of the contribution of the continuous spectrum to h(p, z) in (2.23) is de- 
termined by the support properties of the coupling function f(k); it may contain, in an 
invariant way, restriction to k° positive or negative, or bounds on = — the mass- 
squared of the 9 particle. 

To study these cases, we write (3.2) as 

■»<*> = m~ N - %w e + 557) + s^v^ - p ■ < 3 - 3 > 



14 



where the signs ± in (3.3) correspond to the choice of sign of the energy variable of the 9 
particle, i.e., 

k° = ±^k 2 +m 2 . (3.4) 

The choice of sign determines the direction of propagation of the wave packet of the 9- 
particle in time. As pointed out by Stueckelberg (and Feynman) 9 , the negative energy 
particle has a correspondence with the antiparticle; the CPT conjugate of the wave packet 
is that of a particle with positive energy and the opposite sign of the momentum and the 
charge, i.e., the antiparticle. Applying this definition to a particle with spin 37 , one finds 
that it is completely consistent with Dirac's interpretation (the representations of C, P and 
T on the corresponding Dirac spinors are the same as in the usual on-shell theory) 38 . We 
first discuss the case k° > 0. 

We shall take p° > in the following; the 9 and N particles can exchange roles under 
a slight reformulation, and there is therefore a symmetry in the analysis for p° < 0. Then 
since p° > |p|, if the V particle has positive mass-squared (this condition follows from the 
choice of states of the quantum field), it follows from (3.3) that 

w(k) — > +oo for |k| — > oo. (3-5) 

The lower bound of the spectrum is determined, for each p^,mg : as a function of |k|, by 

dw(k) 
d\k\ 

Let us define the k-dependent part of (3.3) as 

7y+ (k) = p°^k 2 + mj - |p||kj cos 0, (3.6) 
where is the angle between p and k. Then, this condition implies 

r ? ;(k) = ^£L-|p|cos^>. (3.7) 



in 



e 



The minimum value is therefore (the second derivative is positive) achieved for cos > 0, 
at 

p° 2 k 2 = p 2 (k 2 + m^)cos 2 </>, (3.8) 



yielding 



|k| = ; mg, (3.9) 

p° 2 — p 2 COS 2 



so that 

2 _ m d p° 



k 2 + ml = 



\J p° 2 — p 2 cos 2 (J) 



15 



It then follows that 

Mk)) m . n = m*vV 2 -P 2 cos2 0, (3.10) 
which is further minimized by taking (p = 0: 

M k ))„™ = WV- (3-11) 

We note that for cos <f> < 0, the minimum value is at p°mg, always larger than the value 
given by (3.11) for p ^ 0. The spectrum for the two-body N + 9 channel is therefore 



Let us now return to (3.3) for k° < 0. In this case, 

r)-(k) = -p°^Jk 2 + m 2 -p-k (3.13) 

goes to — oo for |k| — > oo, and the spectrum is unbounded from below. The upper bound 
for fixed mg occurs at a value of |k| for which (the second derivative is negative) 

v '_ (k) = ~f |k| = - |p| cos <j>, (3.14) 



yielding, for cos < 0, 



so that 



\A 2 + rrig 



II I P COS ^ (O 1K\ 

|k| = — me, (3.15) 



\fp° 2 ~ P 



2 cos 2 



(l-(^)) raai = -me\Jp° -p 2 cos 2 (/>, (3.16) 
which is further maximized by taking <fi = n: 



(v-(k)) max = -m e ^ (3.17) 
The spectrum for the two body N + 9 channel in the case k < is then 



«,_(*) =(-oo, J^-^(± + J-)- m °^ ). (3.18) 

\ 2M N 2 K M M N > M N ) K 1 

The upper bound of W-(k) is less than the lower bound of w + (k), i.e., 

(w+(k)) min - (w-(k)) max = 2m — — — , (3.19) 

M N 

16 



and there is a gap in the spectrum for every mg 7^ 0. For the case in which mg can reach 
zero, the continuous spectrum, including both k° > and k° < parts, covers the real 
line (—00, 00). 

Let us consider the question, for k° > 0, of whether the "discrete" spectral point 
p 2 /2My is embedded in the continuum. The condition that this point lie above the lower 
bound of w + (k) is that 



rn 



2 



+ i } ~ > T ( iv^ " iv^ } - (3 - 20) 



Although conservation of energy and momentum imply that for a decay process, my > 
rriN + Tag, the dimensional scale fac tors My, Mn, Mg may satisfy this inequality as 
well; when the kinematical masses \f— p 2 , \/—{p — k) 2 and \J —k 2 take on these values, 
the classical squared proper time intervals of the corresponding free world line —dx^dx^ 
become equal to dr 2 ; they may also be taken as the Galilean limit of the kinematical 
masses 12 . In this case, jj^ — is positive. The right hand side of (3.20) is therefore, in 
this case, negative. The minimum value of the left hand side, at 

Mg 



is 



M » >0, (3.22) 



2 M N (M N + Mg) 

and hence, in this case the discrete spectral point is embedded. It is not, however, necessary 
to make the choice My > Mn; a choice of parameters can be made for which the discrete 
state lies below the minimum of w + (k). This would correspond to the isolated "bound 
state" case of the Lee model, for which some studies in mass renormalization have been 
made 39 . 

As for the two body problem discussed in Eq. (1.9) and the following paragraph, the 
kinetic terms of the two body final state can be expressed in terms of the total momentum 
and relative momentum 

= ^- F) + F = |) f (3.23) 

and 

rel M M v ' 

as 

W (k) = — + ?^L, (3.25) 
v 1 2M 2m' V ; 

where M = M N + Mg, m = M N M e /(M N + Mg). 

For the scattering system corresponding to N + 9^N + 9, the Lippmann-Schwinger 
equations provide an effective potential for transitions in p re i, and the problem may be 
treated as discussed in Section 1. This will be done in detail elsewhere. 

For the positive energy contributions, the process V — > + 9 is concerned only with 
particles in the usual sense, and may be represented as a process in space- time as in fig.l. 
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[fig- 1] 

For k < 0, however, the space-time picture appears as in fig. 2, corresponding to 
Stueckelberg's representation of pair annihilation as part of the process V — > N + 9. 

[fig- 2] 

The corresponding physical process, seen in the laboratory, contains the CPT conjugate 
of the 9 external leg, and appears as in fig. 3, corresponding to a real V + 9 annihilation 
to yield the (virtual) N. 

[fig- 3] 

The inverse process, N — > V + 9 corresponds to the physical decay of the N particle. For 
the case n — > p + n~ discussed by Henley and Thirring 27 , this inverse process corresponds 
to p — > n + Tr + . The difficulty that they cite, that in the on-shell relativistic quantum field 
theories both of these processes occur simultaneously in the interaction term of an analog 
relativistic Lee model, does not occur in the model presented here, since the two processes, 
although described in the context of the same theory, occur in two different parts of the 
spectrum. 

For the case k° < 0, the condition that the discrete spectral point is embedded in the 
continuum is 

If My > Mat, the right hand side is negative, and hence the discrete point is an isolated 
bound state. However, in the physical process N — > V + 9, the scale parameters may 
be chosen (as in the argument given for the process V — > N + 9) to satisfy Mn > My, 
in which case the right hand side of (3.23) is positive. The left hand side can extend to 
zero, and hence for this choice, the discrete state is embedded in the continuum as well. 
Since, as we have seen, the inclusion of both positive and negative branches of k° results 
in a spectrum that covers the whole real line (when mg reaches zero), the discrete spectral 
point would then be embedded for any My, Mjq- 

4. The Friedrichs Form and Conclusions 

Since the relativistic Lee model we have considered above has sectors, it can be iso- 
morphically mapped to a spectral model. We write the generator of the motion in the 
direct integral space J» H p over the four- momenta of the V particle as (for each p^) 

K p = Mvip)Po + K p + K I , (4.1) 

where Ki has matrix elements only between the continuum and the discrete state cf> p 
(here, P = <j> p (l>p) with eigenvalue My(p) corresponding to p 2 /2M V: and K p has absolute 
continuous spectrum (with, perhaps, finite or semi-infinite gaps). This form corresponds 
to the canonical decomposition 

K p = P K p P + PK p P + P K p P + PK p P , (4.2) 

where P is the projection into the continuous spectrum (at each point p). 

The vertex function f(k) of the field theoretical model can be covariantly restricted 
to support in the upper light cone, k° > 0, and the spectrum (with bounded mg) will be 
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therefore similar to that of the usual non-relativistic Friedrichs model, i.e., for which the 
unperturbed evolution operator has semi-infinite absolutely continuous spectrum and an 
embedded discrete state with an eigenvalue corresponding to the mass of the initial state 
of the decay problem (note a similar phenomenon in structure of the Feynman propagators 
studied in ref. 40). 

As we have seen, if k can be positive or negative, the continuous spectrum can reach 
±00; if me can reach the value zero, the continuous spectrum will have no gaps, and the 
spectral model is very similar to that of Pietenpol 41 ; if me is bounded away from zero, 
there will be a gap. Furthermore, the discrete state may occur embedded in either branch 
of the continuum or may occur within the gap. There are therefore a large variety of 
covariant models which may be extracted from the framework given here. 

Denoting, for each p, the second term on the right of (4.1) by 

K p = J w(k)\k)(k\d 4 k, (4.3) 

a calculation of the type given above, from (1.28) to (1.35), results directly in the relation 
(2.23) for the reduced resolvent at each p. 

Somwhat more abstractly, introducing a continuous label A, we may denote the con- 
tinuum generalized eigenstates by {|A)}, so that (4.3) takes the form 



K p = J w(X)\X)(X\d\. 



The reduced resolvent then takes the generalized Friedrichs form 

h(z,p) = z-Mv(p)- [dX K^lMl! . (4.4) 

J z — w(X) 

The complex poles which dominate the decay law (in "proper time"), as discussed at the 
end of Section 2, can then be investigated in a way similar to the discussion of the nonrel- 
ativistic form in Section 1. The construction of the generalized states 4 of the relativistic 
Lee-Friedrichs model can then be carried out. 

We remark here that the generalized eigenstate in this construction corresponds to 
a shift is the position of the unperturbed eigenvalue of K p from A4v(p) to the complex 
zero of h(p, z) in the second Riemann sheet. This can be interpreted as the acquisition 
of a complex part to the total energy- momentum of the system p^. Let us label the pole 
position by (we neglect the dependence of M', 7 on p here) 

M' 2 7 

C=-— i-, (4-5 

S 2M V 2' V ; 

so that, as for the nonrelativistic case, the decay law goes as e~ JT . Then, in the center of 
momentum frame, 

p 2 E 2 



2M V 2M V 



C, 
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or, for 7 < M', 



(4.6) 



For M' = 



My + 5My ( we assume 




for the unperturbed system), 




(4.7) 



where we have neglected a term of order ^(SMv /My)- The result (4.7) illustrates how, in 
the rest frame of the initial particle, the mass change and decay width generated by the 
interaction in this covariant model can emerge as a complex energy. 
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Figure Captions 

fig. 1. The process N -> V + 6(k° > 0) 
fig. 2. The process N -^V + 9(k° < 0) 
fig. 3. The process N + 9(k° > 0) -> ^(virtual) 
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